Reaching the continuum limit in lattice gauge theory - without a computer. 
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The scaling slope of the anti-symmetric mass gap M of compact U(I)2+i lattice gauge theory 
is obtained analytically in the Hamiltonian formalism using the plaquette expansion. Based on the 
first four moments of the Hamiltonian with respect to a one-plaquette mean field state the results 
demonstrate clear scaling of M at and beyond the transition from strong to weak coupling. The 
scaling parameters determined agree well with the range of numerical determinations available. 



PACS number: 11.15.Ha 

With the vast majority of calculations carried out in 
lattice gauge theory being of an intensely numerical na- 
ture, non-perturbativc analytic results are of considerable 
interest. If we look to complex theories such as QCD, 
the most basic quantities one wishes to calculate are the 
states of the mass spectrum of the theory. In lattice 
gauge theory, such non-perturbative calculations are car- 
ried out numerically in the path integral representation 
to some level of approximation by evaluating the correla- 
tion functions by Monte-Carlo simulation. Typically, one 
must work with large lattices and spacings which are fine 
enough to detect the continuum behaviour of the the- 
ory in question. The computational resources employed 
in this approach are enormous, essentially reflecting the 
seemingly intractable many-body problem at the core of 
non-perturbative quantum field theory. We report here 
on a novel approach to studying lattice gauge theories, 
and as an indication of its potential apply the method 
to compact U(l) lattice gauge theory in 2+1 dimensions. 
This theory, although obviously simpler than QCD it- 
self, nevertheless shares many of the features of the QCD 
many-body problem on the lattice, namely similar nu- 
merical difficulties in calculating masses in the contin- 
uum limit. The scaling behaviour for the mass gap, M, 
is expected to be: 
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where fco and fci are constants, a is the lattice spacing 
and f3 = 1 / g 2 where g is the dimensionless coupling con- 
stant. As with most other interesting lattice gauge the- 
ories, it has long been a maxim that even U(l)2+i lat- 
tice gauge theory must be studied by large scale numer- 
ical computation in order to determine the lowest mass 
state of the theory. Indeed, although studied by a range 
of different methods, the scaling properties of this the- 
ory are still not known precisely: fco ~ 4.1 — 6.345 and 
k\ ~ 4.369 — 6.27 (see M for summary). In this sense the 
results we present are quite spectacular - an analytic de- 
termination of the mass-gap, which clearly demonstrates 
the scaling behaviour expected in the continuum limit 
with fc = 4.42 and k% = 5.41. 



Throughout this paper we will work in the Hamilto- 
nian formalism. When one speaks of "lattice gauge the- 
ory" it is usually meant the Lagrangian formalism where 
one essentially converts the theory to a statistical model 
in Euclidean space into order to facilitate the use of pow- 
erful Monte-Carlo methods. In contrast, there is a dearth 
of equivalently efficient methodology in the Hamiltonian 
formalism of lattice gauge theory. This fact has been re- 
sponsible for the domination of action-based lattice gauge 
theory over the Hamiltonian approach. While there ex- 
ist numerical and Monte-Carlo methods applicable to the 
Hamiltonian lattice formalism, these have been of limited 
use, and certainly the results are not of the same quality 
as in the action-based approach. However, there is a cer- 
tain appeal of the Hamiltonian formalism which tugs at 
physicists' intuition. There are also technical reasons for 
considering the Hamiltonian approach - one often works 
directly on infinite lattices, fermions are simpler to incor- 
porate (no fermion determinant), the lattices are 3D as 
opposed to 4D (time is continuous), and finite temper- 
ature physics may actually be more natural in this for- 
malism [pi . The lack of powerful and systematic method- 
ology, such as is available in the Lagrangian formalism, 
has prevented these features from being exploited or in- 
vestigated - hence, the importance of new methodologies. 
Here we use a novel approach - the plaquette expansion 
- based on a large volume expansion of Lanczos tridiag- 
onalization. 

The Kogut-Susskind Hamiltonian for U(l)2+i on a lat- 
tice of N plaquettes is given by B : 
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The strong-coupling limit is defined by g — * oo and the 
weak-coupling limit by g — > 0. For the U(l) model, Ei is 
the electric field operator associated with lattice link I. 
In the representation in which Ei is diagonal, the eigen- 
value Ei takes on integer values 0,± 1, ± 2,... on each 
link and the link operators Ui and (jj raise or lower the 
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value of Ei on link I by unity respectively due to the 
commutation relations: 



[E u Ui] = Ui 



[Ei,U, 
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The operator U p acts on the links around the smallest 
closed (Wilson) loop or square on the lattice known com- 
monly as a plaquette p: 



U p = u&ul&l 



(4) 



The extensive nature of the Hamiltonian Equation (g) 
lends itself to the application of the plaquette expansion. 
Indeed this method was originally devised to study such 
local lattice models ||. Beginning with the well known 
Lanczos recurrence for the iterative transformation of the 
Hamiltonian into tri-diagonal form, a trial state \ip±) with 
the desired symmetries of the state of interest |\&o) (the 
ground state or vacuum in this case) and non-zero over- 
lap (i.e. (-01 l^o) ^ 0) is chosen and a basis {|^ n )} is 
constructed according to 
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where a n = {ip n \H\ip n ) and j3 n = (tp n +i\H\ip n } are the 
matrix elements of the Hamiltonian in tri-diagonal form. 

The matrix elements, a n and /3 n , can be written in 
terms of Hamiltonian moments (H n ) = (ip 1 \H n \ipi). As 
the connected part of the Hamiltonian moment is pro- 
portional to the volume of the system (H n ) c = e n N, 
one may utilize the cumulant transformation ||, 



(H n ) 



E 

p=0 



^ (HP+ 1 ) c {H n - 1 -P) 



(0) 



and rewrite the a n and (3 n in terms of the connected vac- 
uum coefficients e n . Extensivity of the problem leads to 
the following plaquette expansions in 1/N [||: 
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where s = n/N. In the bulk limit (n, N — > oo), keeping 
s fixed, it was shown in |t]] that one may perform the 
exact diagonalization of the Lanczos tri-diagonal matrix 
for the vacuum energy density £q: 



£ = inf (a(s) - 2(3(s)) 



(9) 



For low order truncations of the series (employing mo- 
ments up to (H r )) one can minimize (0) with respect to 
s analytically and obtain the closed form expressions for 
the ground state energy, So(r): 

£o(l) = ei, 
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Note the first approximant £o(l) is just the variational 
result. The expressions for £o(3) and £o(4) clearly illus- 
trate the non-perturbative nature of the plaquette expan- 
sion as an effective summation over clusters of increasing 
complexity. 

To determine expressions for excited states we note 
that lowest state of a given sector, S, in the spectrum 
will in general have an energy: 
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where MX is the mass gap. For excited states the 
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Hamiltonian moments have the form (H r 

(s) 

m n . Thus, we can then define analogous "cumulants" 

el 5 ^ for this sector of Hilbert space as = e n +rrin /N, 
(s) 

where we call m n the mass gap connected moments. 

{ s) 

One can derive expressions for approximants, Mq (r), 
to the mass gap, in a similar fashion to those for the vac- 
uum energy density. For example, using up to 4th order 
moments we have [pi : 
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where the vacuum moment functions, F ni are given by 
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with A = 3 e 3 2 — 2 e 2 e 4 and = 8 e 2 e^ 2 e 4 — 6 — 
3e 2 2 e 4 2 . 

It is appropriate to comment on relevant methods, 
which by and large employ trial states encompassing as 
much of the expected continuum physics as is viable. Of- 
ten, such calculations have employed the simplest gauge 
invariant trial state — the strong-coupling vacuum |0) 
(defined as the state satisfying Ei\0) = for every link). 
This is the perturbative starting point for series calcula- 
tions, which are then extrapolated to weak coupling. Al- 
though the calculation of moments with respect to this 
state can be carried out to relatively high order, typically, 
one finds that this state is simply inadequate to explore 
the non-perturbative weak-coupling regime of the the- 
ory. In the case of the plaquette expansion, a window of 



scaling in terms of the expansion order was evident, but 
higher order results became problematic |j) . 

More recently, trial states which are more complex and 
appropriate to weak-coupling have been investigated by 
a number of authors JTo[-p~5| in varying contexts. There is 
a clear tradeoff between physical complexity of the trial 
state, and algebraic complexity of the moment calcula- 
tion. Given this, we utilized a one plaquette exponen- 
tial trial state, that has the correct weak coupling be- 
haviour for the vacuum energy density - the so-called 
mean-field state - defined as: \ipi) = exp{S f /2}|0), where 
S = A E P {^p+^p} an d A will eventually be determined 
variationally from (H) . It is convenient to work in the so- 
called U representation where the trial state is a function 



of Up variables: 
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As the state is constructed from plaquette variables, U p , 
it is automatically gauge invariant. One expects the true 
ground state of the system to have contributions from 
all possible Wilson loops - an important feature of the 
plaquette expansion method is that by its very nature 
larger sized Wilson loops (larger clusters) are systemati- 
cally introduced. 

In D = 2 + 1 the Hamiltonian moments can be written 
in terms of integrals over plaquette variables, for example 
the vacuum moments are given schematically by: 
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For the group U(l), the plaquette variables may be 
rewritten in terms of plaquette phase angles U p = 
exp(i<p p ). The calculation of Hamiltonian moments for 
use in the plaquette expansion utilizing the one-plaquette 
state in 2 + 1 dimensions comes down to evaluating inte- 
grals of the following form, 
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where Jo(2A) is the zeroth-order modified Bessel function 
of the first kind. The integrand functions, G n [4> p \, natu- 
rally become more complex as one increases the moment 
order, however, all integrations can be carried out ana- 
lytically. All expressions are analytic in A and (3, and 
once (H) is minimized we implicity determine A(/3). 

It is well known |l2|,[l6| that the variational calcula- 
tion (H) applying the one-plaquette state in D = 2 + 1 
dimensions involves no correlations between plaquettes 
and thus one is simply solving a Mathieu-type problem. 
In this context, the benefit of choosing the one-plaquette 
state as the trial state rather than the strong-coupling 
vacuum becomes more apparent on investigation of the 



strong (/?— ►O) and weak-coupling (/3— >oo) limits of (H). 
For the strong-coupling limit, one finds for the energy 
density: e (l) = (H)/N ~ j3 - P 3 /A + 0(/3 7 ) , which 
agrees with strong-coupling perturbation theory [ jl7| . 
Similarly, for the weak-coupling limit, one finds: £o(l) ~ 
A-A 2 /(8(3) + 0(l//3 2 ) , with A = 1. This agrees exactly 
with weak-coupling perturbation theory for the Mathicu 
problem to the order shown. Exact U(l)2+i has a similar 
weak-coupling expansion |l8| , |l9fl with A = 0.9581. With 
the inclusion of larger sized Wilson loops obviously the 
U(l) calculations, as they should, will deviate from the 
Mathieu results. As this trial state encompasses the ac- 
tual behaviour in both limits, one has confidence in using 
the plaquette expansion to systematically obtain results 
throughout the range of (3. 
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FIG. 1. Mass gap results applying one (minimized numeri- 
cally), three and four one-plaquette Hamiltonian moments. 
Scaling is evident for the four order moments curve from 
(3 = 0.7 to 1.5. 

Using a diagrammatic method (details to be presented 
elsewhere) we were able to derive the first four one- 
plaquette connected moments for both the vacuum and 
anti-symmetric state. These moments immediately give 
analytic approximates for vacuum energy and the mass 
gap. 

Since the expected scaling behaviour for compact 
U(l)2+i is of the form (Q), the one-plaquette state results 
for the plaquette expansion to first order are presented 
as a In (M 2 //3) versus (3 plot in Figure [|. The results us- 
ing one, three and four moments appear to be converging 
towards the scaling slope. Clearly scaling is evident for 
inverse-coupling values (3 = 0.7 to 1.5, passing the transi- 
tion point at (3 = 1. The scaling form we find is given by 
M 2 — (3 exp (—4.34/3 + 5.34), and agrees well with other 
estimates (as summarized in 0). 

It should be noted that the minimization of (H) is the 
only strictly numerical procedure (small as it is) in our 
calculation. A totally analytic estimation (keeping in line 
with the title of this paper) is possible when one writes 
the function A(/3) which minimizes (H) as a series expan- 



sion: A = /3 2 /2— /3 6 /16+. . . . Keeping the first five terms 
of the expansion (for which the resulting values have con- 
verged), scaling is observed from j3 — 0.8 to 1.2 (Fig- 
ure |). The series breaks down at larger (3 values. The 
scaling form is given by M 2 — (3 exp (—4.42/3 + 5.41), and 
again agrees well with the numerical estimates available. 

In summary, applying the plaquette expansion we were 
able to identify scaling of the anti-symmetric (or pho- 
ton) mass gap in the weak-coupling regime of compact 
U(l)2+i lattice gauge theory via a completely analytic 
calculation utilizing only four one-plaquette state Hamil- 
tonian moments. Furthermore, although we have pre- 
sented results for only the first four moments (for which 
the calculation is entirely analytic) it should in principle 
be possible to derive higher moments with respect to this 
state. Indeed wc have preliminary results for the vacuum 
energy density to sixth order Q giving high accuracy over 
a large range of couplings. 
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FIG. 2. Convergence of the complete analytic mass gap 
results for one (variational), three and four one-plaquette 
Hamiltonian moments to the continuum limit, using the series 
form for A(/3). 

These results may have important implications for lat- 
tice gauge theory - the plaquette expansion method itself 
can be applied to any extensive many-body problem. The 
only existing application of the method to QCD3 + i (in- 
cluding massless quarks) with respect to the strong cou- 
pling state proved somewhat disappointing pQ ] in that 
scaling of the nucleon/meson mass ratio was not evident. 
However, the results presented here indicate that the one- 
plaquette state may be of sufficient quality to capture the 
essence of the continuum physics in lattice gauge theory 
at relatively low orders. An exciting possibility is that the 
continuum limit might finally be observed in the Hamil- 
tonian formalism for QCD using this method. 

To finish, we propose a semi-analytic program to study 
the 3+1 lattice gauge theories in the Hamiltonian formal- 
ism. The main obstacle to overcome is the calculation 
of the moments in 3+1 dimensions. The appeal of 2+1 



systems as test models is that the transformation from 
link variables to plaquette variables is trivial and makes 
the integrations tractable - hence the analytic work re- 
ported in this paper. In 3+1 dimensions it is well known 
that the transformation cannot be carried out in closed 
form due to the appearance of Biachi identities. How- 
ever, the calculation of moments in 3+1 dimensions by 
Monte-Carlo methods is actually a relatively small scale 
numerical exercise. Because we are dealing with a cluster 
expansion (we need only connected moments) the lattices 
required fairly modest in extent, and are by definition 
only three dimensional. The integrands become compli- 
cated as larger correlations are included, which means 
that the statistics must be very good, however, prelim- 
inary calculations have demonstrated that the required 
precision is possible to achieve with quite modest com- 
puting resources 0. 
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